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Abstract 

We show that the gravitational modification of the phase of a neutron beam (the 
COW experiment) has a classical origin, being due to the time delay which classical 
particles experience in traversing a background gravitational field. Similarly, we show 
that classical light waves also undergo a phase shift in traversing a gravitational field. 
We show that the COW experiment respects the equivalence principle even in the 
presence of quantum mechanics. 

In a landmark series of experiments Q, || Colella, Overhauser and Werner (COW) and 
their subsequent collaborators (see e.g. Refs. |J, Q for overviews) detected the modification 
of the phase of a neutron beam as it traverses the earth's gravitational field, to thus realize 
the first experiment which involved both quantum mechanics and gravity. A typical generic 
experimental set up is shown in the schematic Fig. (1) in which a neutron beam from a reactor 
is split by Bragg or Laue scattering at point A into a horizontal beam AB and a vertical 
beam AC (we take the Bragg angle to be 45° for simplicity and illustrative convenience in 
the following), with the subsequent scatterings at B and C then producing beams which 
Bragg scatter again at D, after which they are then detected. If the neutrons arrive at A 
with velocity Vq (typically of order 2 x 10 5 cm sec -1 ) and ABCD is a square of side H, then 
the phase difference 4>cow — 4>acd — 4>abd is given by —mgH 2 /hv to lowest order in the 
acceleration g due to gravity |] , and is actually observable despite the weakness of gravity, 
since even though / p-df only differs by the very small amount i7i(vcd~vab)H = —mgH 2 /vo 
between the CD and AB paths, nonetheless this quantity is not small compared to Planck's 
constant, to thus give an observable fringe shift even for H as small as a few centimeters. 

The detected COW phase is extremely intriguing for two reasons. First, it shows that it 
is possible to distinguish between different paths which have common end points, with the 
explicit global ordering in which the horizontal and vertical sections are traversed leading to 
observable consequences. And second, it yields an answer which explicitly depends on the 
mass of the neutron even while the classical neutron trajectories (viz. the ones explicitly 
followed by the centers of the wave packets of the quantum mechanical neutron beam) of 
course do not. The COW result thus invites consideration of whether the detected ordering 
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is possibly a topological effect typical of quantum mechanics, and of whether quantum me- 
chanics actually respects the equivalence principle.^ As we shall see, the ordering effect is 
in fact already present in the motion of classical particles in gravitational fields, and even in 
the propagation of classical waves in the same background, with this latter feature enabling 
us to establish below that the mass dependence of the neutron beam COW phase is purely 
kinematic with the equivalence principle then not being affected. 

To address these issues specifically we have found it convenient to carefully follow the 
neutron as it traverses the interferometer, to find that the two beams do not in fact arrive at 
the same point D or even at the same time, with this spatial offset and time delay not only 
producing the interference effect, but also being present in the underlying classical theory. 
Quantum mechanics thus does not cause the time delay, rather it only serves to make it 
observable. Since gravity is a relativistic theory we shall need to introduce curvature (which 
we do below), but we have found it more instructive to consider the non-relativistic limit first. 
Since we can treat the neutron beams as rays, their motions round the ABCD loop can be 
treated purely classically between the various scatterings. Moreover, the various scatterings 
themselves ai A, B, C and D introduce no additional phases, are energy conserving, and 
give angles of reflection equal to the angles of incidence.^ Thus the entire motion of the 
neutron is the same as that of a classical macroscopic particle (the neutron spin plays no 
explicit role in the COW experiment, so we treat the neutrons as spinless) which undergoes 
classical mirror reflections. 

The neutron which goes up vertically from A to C arrives at C with a velocity (0, vq — 
gH/vo). The neutron AC travel time is t(AC) and the classical action S — f(p ■ df — E$dt) 
(Eq = mvl/2) undergoes a change S(AC) where 

t(AC) = H/v + gH 2 /2vl , S(AC) = mv H - mgH 2 /2v - E Q t(AC) (1) 

On scattering at C the neutron is then reflected so that it starts off toward D with a velocity 
{vq — gH/v ,0). On its flight it dips slightly to arrive at the next scattering surface at the 
point Di with coordinates (H — gH 2 /2v^, H — gH 2 /2vq), so that there is a change in the end 
point of the motion which is first order in g and thus relevant to our discussion. At D\ the 
neutron has a velocity (t> — gH/vo, —gH/v ), with the CD\ segment taking a time t(CD\) 

1 Possible implications of the mass dependence of the COW phase on another related issue, viz. quantum 
mechanics and the action and reaction principle, are explored in Ref . |J . 

2 This is not completely obvious since if the neutron slows down a little or is deflected a little, then at 
each subsequent scattering it is not quite incident with a wavelength and angle which satisfy a Bragg peak 
condition. However, consider a crystal for which Bragg scattering would occur for the pair of vectors ki and 
kf, viz. equal length vectors which obey fcj — kf = G where G is a reciprocal lattice vector which imparts a 
momentum transfer perpendicular to the scattering surface. Now instead allow ki to impinge upon a crystal 
with additional momentum Afcj with the outgoing kf beam then emerging with additional momentum Akf. 
Since each such scattering is coplanar the most general such momentum modification can conveniently be 
written as Afcj = onki + afkf, Akf = (3iki + with energy conservation yielding fii = otf, (if = on since 

there is no energy transfer to the crystal. Consequently, the momentum transfer at a scattering is now given 
by ki + Aki — kf — Akf = (1 + aj — (Xf)G, and is thus still perpendicular to the scattering plane, but now no 
longer right on a diffraction maximum. The angle of reflection is thus still equal to the angle of incidence, 
but with a scattering intensity which is slightly reduced with respect to the Ak = case. (The only way 
that we would actually be able to shift into a new Bragg peak condition would be to have a,; — a/ = 0, a 
condition which only occurs when each Afc is parallel to the scattering surface.) 
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and contributing an amount S{CDi) to the classical action where 

t(CDi) = H/v + gH 2 /2vl , S(CT>i) = mw if - 3mgH 2 /2v - E^CD^ (2) 

The neutron which starts horizontally from A arrives not at B but at the point B\ with 
coordinates (H — gH 2 /2v 2 , —qH 2 /2vq) and with a velocity (vq, —gH/vo). The AB\ segment 
takes a time t(ABi) and the action changes by S(ABi) where 

t(ABi) = ff/uo - gH 2 /2vl , = m^tf - mgH 2 /2v - E^AB^ (3) 

After scattering at i?i the neutron sets off toward D with velocity (—gH/vo,vo) and arrives 
not at D or £>i but rather at the point D 2 with coordinates (H — 3gH 2 /2v 2 : H — 3gH 2 /2vq), 
and reaches there with velocity (—gH/vo,vo — gH/vo). The B\D 2 segment takes a time 
t(BiD 2 ) and the action changes by S(BiD 2 ) where 

t(B 1 D 2 )=H/v -gH 2 /2v 3 , S(B 1 D 2 ) ^ mv H - 3mgH 2 /2v - Eo^D^ (4) 

As regards the neutron's path around the loop, we see that the small vertical gH 2 /2v 2 dip 
during each of the two horizontal legs causes each neutron to travel a distance gH 2 /2v 2 less in 
the horizontal than it would have done in the absence of gravity, to thus provide a first order 
in g modification to / p-df in each of these legs, even while these same vertical dips themselves 
only contribute to the action in second order. However, for the two horizontal sections, each 
leg is shortened by the same amount in the horizontal, so that the difference in / p ■ df 
between the CD 1 and AB 1 legs still takes the value m(v C D — vab)H = —mgH 2 /v quoted 
earlier. As regards the two vertical legs, we note that even though the AC leg is completely 
in the vertical, since the neutron starts the B\D 2 leg with a small horizontal velocity, during 
this leg the neutron changes its horizontal coordinate by an amount gH 2 /v 2 , thereby causing 
it to reach D 2 after having also traveled a distance gH 2 /v 2 less in the vertical than it would 
travel in the AC leg. Consequently, there is both a spatial offset (gH 2 /v 2 , gH 2 /vl) between 
D 1 and D 2 , and a time delay t(AC) + i(CDi) - t(ABi) - *(-Bi-D 2 ) = 2 gH 2 /v% between the 
arrival of the two beams, with the integral j p ■ df thus not taking the same value in each 
of the two vertical legs. However, our calculation shows that all the these modifications 
actually compensate in the overall loop with there being no difference in J p-df between the 
ACDi and AB\D 2 paths. However, since there is a time delay, there is still a net change 
in the action, given by S(AC) + S(CL>i) - S(ABi) - S^B^) = -mgH 2 /v . However, 
we cannot identify this quantity with the COW phase hA<f)cow since the beams have not 
interfered due to the spatial offset between D\ and D 2 . 

Before discussing the issue of this spatial offset, it is instructive to ask where the neutrons 
would have met had there been no third crystal at D to get in the way. Explicit calculation 
shows that they would in fact have met at the asymmetric point D 3 with coordinates (if — 
3gH 2 /2v 2 , H-gH 2 /2vl) with the CD 3 and £iD 3 segments taking times t(CD 3 ) and t( J Bi J D 3 ) 
while yielding action changes S(CD 3 ) and S(BiD 3 ), where 

t(CD 3 ) = H/v - gH 2 /2vl , S(CD 3 ) = mv H - 5mgH 2 /2v - E t(CD 3 ) 
t(B 1 D 3 )=H/v + gH 2 /2v 3 , S^-Bi-t^) = mv H — mgH 2 /2v — E Q t(BiD 3 ) (5) 

The neutrons would thus meet at D 3 without any time delay, and with S(AC) + S(CD 3 ) — 
S(AB 1 ) — S(B 1 D 3 ) = —2mgH 2 /v . We thus see that for purely classical particles reflecting 
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off mirrors at B\ and C the quantity / p ■ df evaluates differently for the two paths ACD3 
and AB1D3. This is a thus a global, path dependent effect in purely classical mechanics 
in a background classical gravitational field which is completely independent of quantum 
mechanics.^ However, since the classical action is not observable in classical mechanics, it 
is only in the presence of quantum mechanics that phase differences become observable. (In 
classical mechanics what is observable is that the neutrons meet at D 3 rather than on the 
AD axis.) 

Returning now to the COW experiment itself, in order to understand the implications 
of the time and spatial offsets between D\ and D2, it is instructive to consider the Young 
double slit experiment with purely classical light. As shown in Fig. (2), light from a source 
S goes through slits Q and R to form an interference pattern at points such as P, with 
the distance Ax = QT representing the difference in path length between the two beams. 
Given this path difference, the phase difference between the two beams is usually identified 
as kAx, from which an interference pattern is then readily calculated. However, because of 
this path difference, the SQP ray takes the extra time At = Ax/c to get to P, to thus give 
a net change in the phase of the SQP beam of kAx — uAt which actually vanishes for light 
rays. The relative phase of the two light rays in the double slit experiment thus does not 
change at all as the two beams traverse the interferometer. However, because of the time 
delay, the SRP beam actually interferes with an SQP beam which had left the source a 
time At earlier. Thus if the source is coherent over these time scales, the SQP beam carries 
an additional +uAt phase from the very outset. This phase then cancels the —uAt phase it 
acquires during the propagation to P (a cancellation which clearly also occurs for quantum 
mechanical matter waves moving with velocities less than the velocity of light), leaving just 
kAx as the final observable phase difference, a quantity which is non-zero only if there is in 
fact a time delay. We thus see that the double slit device itself actually produces no phase 
change for light. Rather, the choice of point P on the screen is a choice which selects which 
time delays at the source are relevant at each P, with the general interference pattern thus 
not only always involving the time delay at the source, but also in fact always requiring one. 

With this in mind, we now see that we also need to monitor the time delay of the neutron 
in the COW experiment, However, since the total energy of the neutron does not change as 
it goes through the interferometer, the time delay contribution will still drop out of the final 
phase shift expression (explicitly but not implicitly). However, for the COW experiment we 
noted above that as well as a time delay between the ACD\ and AB1D2 paths, there was 
also a spatial offset. Consequently, the AR>iD 2 path interferes not with the ACD\ path, but 

3 Since the momentum p is the vector derivative VS 1 of the stationary classical action, for any trajectory 
which obeys the equations of motion the integral J p ■ df then only depends on its endpoints (in a manner 
which in general depends on the explicit VS* direction in which each endpoint is approached.) Consequently, 
it is at first a little surprising that J p ■ df could evaluate differently for the ACD3 and AB1D3 paths, 
since these paths do share common endpoints. However, neither one of these paths is the one which is the 
stationary one for the motion which starts at A and finishes at D3 (viz. a parabola close to the straight 
line ^4.D3); and thus even while the ACD3 path for instance is composed of the two subpaths AC and CD3 
each one of which is stationary between its own endpoints, nonetheless the sum of the two subpaths is not 
stationary between the overall A and D% endpoints, with p then not being expressible as a total derivative for 
the whole ACD3 (or AB\D%) motion. In passing, we note that the fact that the sum of two stationary paths 
is not necessarily stationary between the overall endpoints is also the origin (see e.g. [||) of the existence of 
the many paths of the Feynman path integral description of quantum mechanics. 
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rather with the indicated offset A\C\D 2 path, a very close by path which in fact is found to 
lie a distance gH 2 /vq vertically below AB, an offset distance which is within the resolution 
of the beam. The evaluation of the phase shift is then exactly as before with S(AiCiD 2 ) 
taking the exact same value as S(ACDi) to lowest order in g. Now we noted above that all 
the J p-df contributions actually cancel for this particular set of paths. However, because of 
the spatial offset between Dx and D 2 , the AB\D 2 path beam has to travel an extra horizontal 
distance A 2 A = gH 2 /v^ to first get to the interferometer (to therefore provide an analog to 
the distance Ax = QT in the double slit experiment, with A\ and A 2 acting just like the 
pair of slits Q and R). Now in traveling this extra A 2 A distance this beam actually acquires 
yet another time delay t(A 2 A), to therefore impose yet another relative phase condition at 
the source which then identically cancels the associated — E t(A 2 A) change in the action. 
Moreover, in traveling this extra A 2 A the integral J p-df acquires yet one more contribution 
mgH 2 /vQ, and this term then emerges as the only contribution in the entire circuit which 
is not canceled. Consequently, we obtain A<pcow — —mgH 2 /hvo as the final observable net 
COW phase shift. f] 

Turning now to a fully covariant analysis, we need to look at solutions to the Klein- 
Gordon equation (j)^ — (mc/h) 2 <j) = (0 M denotes dcft/dx^) in generic background fields of 
the form dr 2 = B(r)c 2 dt 2 — dr 2 /B(r) — r 2 dVt where B(r) = 1 — 2MG/c 2 r. Firstly we note 
that the non-relativistic reduction of this Klein-Gordon equation is straightforward, with the 
substitution =exp(—imc 2 t /h)ip then yielding 

.. dib h 2 _„ , mc 2 , , . „. mMG 

for slowly moving particles. We thus see that the inertial mass m which is defined via 
the starting Klein-Gordon equation thus also serves as the passive gravitational mass which 
serves to couple massive particles to gravity, so that the particle modes associated with the 
quantization of the Klein-Gordon field thus automatically obey the equivalence principle, 
precisely because of quantum mechanics in fact.f] 

As regards the covariant Klein-Gordon equation, it is convenient to make the substitution 
(f)(x) =exp(iS(x)/h), so that the phase then obeys S^S^ + m 2 c 2 = ihS^.^. In the eikonal or 
ray approximation the ihS^ term can be dropped, so that the phase S(x) is then seen to 

-s fu A n t i/'n i -."1 ficivi ) 



obey the classical Hamilton- Jacobi equation S^Su + m 2 c 2 = 0, an equation whose solution 



4 While we obtain the same answer for A(f>cow as previous authors, nonetheless our derivation seems to 
be somewhat different from the previously published ones. 

5 We note that since the starting covariant Klein-Gordon equation only possesses one intrinsic mass scale, 
it would appear that the equivalence principle has to emerge. However, there is actually a hidden assumption 
in the use of the covariant Klein-Gordon equation (see since it is not the most general equation in curved 
space which can reduce to the flat Klein-Gordon equation in flat space. Rather, the curved space equation 
could also possess additional explicitly curvature dependent terms, terms which would then vanish in the 
flat space limit, but which would modify the particle's coupling to gravity in curved space in a potentially 
equivalence principle violating manner. Now such possible terms (which would of course have to produce 
effects of order less than 10 -11 so as not to spoil the Eotvos result) are simply ignored in standard gravity 
without any apparent justification as far as we can tell, with the standard gravitational phenomenology 
only in fact following in their assumed absence. Thus in passing it is of interest to note that any such 
possible additional curvature dependent terms are not in fact allowed to appear Q in the conformal invariant 
gravitational alternative currently being considered by Mannheim and Kazanas (see e.g. [0). 
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is the stationary classical action between relevant end points. We thus establish that in the 
eikonal approximation the phase of the wave function of a material particle is in fact the 
classical action just as required for the discussion of the COW experiment we gave earlier. 
In the eikonal approximation we can also identify as the momentum p^ = mcdx^/dr, so 
that we can set S(x) = J p^dx**; with the covariant differentiation of the Hamilton- Jacobi 
equation then yielding p^p 1 ' = 0, which we recognize as the massive particle geodesic 
equation. 

In order to actually calculate the geodesies in the gravitational field of the earth it is 
convenient to rewrite the Schwarzschild metric in terms of a Cartesian coordinate system 
x = rsin8cos(f), y = rsinOsirup, z = rcosO — R erected at a point on the surface of the earth. 
With z being normal to the earth's surface, to lowest order in x/R, y/R, z/R, MG/c 2 R 
(where M is the mass of the earth and R its radius) the Schwarzschild metric is then found 
to take the form 

dr 2 = f(z)c 2 dt 2 -dx 2 -dy 2 -dz 2 /f(z) (7) 

where f(z) = 1 — 2MG/c 2 R + 2gz/c 2 and where g denotes MG/R 2 . For this metric the 
non-relativistic geodesies for material particles are given by 'x = 0, y = 0, z = — g, to thus 
enable us to completely justify our earlier non-relativistic calculation.^ 

As regards the purely classical, massless case, on defining <j){x) —exp(iT(x)), we can this 
time identify the eikonal phase derivative T M with the wave number = dx^/dq where q is 
a convenient affine parameter which can be used to measure distances along trajectories. In 
the massless case the Hamilton- Jacobi equation takes the form 

f(z)(k ) 2 - (k 1 ) 2 - (k 2 ) 2 - (k 3 ) 2 /f(z) = (8) 

and again yields the requisite massless particle geodesic equation k^k v = just as in the 
massive case. In the Schwarzschild geometry of interest these massless geodesic equations 
are found to take the form 

A; = cdt/dq = a / f(z) , k 1 = dx/dq = a\ 
k 2 = dy/dq = a 2 , k 3 = dz/dq = (a 2 - (a\ + al)f(z)) l/2 (9) 

where the a« are integration constants. From these equations we recognize the existence of 
the gravitational frequency shift (since k° depends on z), the gravitational time delay (dz/dt 
depends on z), and the gravitational bending of light (dz/dx depends on z if a,\ ^ 0). All 

6 Since the metric of Eq. (|?|) yields a uniform gravitational acceleration, it must therefore be flat, and 
indeed it can readily be transformed into a flat Cartesian metric dr 2 = c 2 dt' 2 — dx 2 — dy 2 — dz' 2 via the 
coordinate transformation || ct' — c 2 sinh(gt/c)f 1 / 2 (z)/g 1 z' — c 2 (cosh(gt/c)f 1 ^ 2 (z) — f 1 ^ 2 (0))/g under 
which the small (z,t) region transforms into the small (z',t') region. However, under such a transformation 
the time delay in the COW experiment would not be removed (the time delay between the two beams 
at the source is a timelike vector which remains timelike under any coordinate transformation), with the 
fact of a time delay thus being a covariant indicator for the COW effect. Similarly, since the action is the 
covariant scalar J p^dx 11 , its value in any trajectory must remain unchanged under this transformation to 
flat Cartesian coordinates. Now while there is no explicit reference to g in flat coordinates, nonetheless, in 
an accelerated coordinate frame in flat space there is still an implicit dependence on g, since points such as 
D will accelerate upward with acceleration g causing the beam from C to arrive at the g dependent D\ , to 
thus expressly show the need for first order changes in the end points of the interferometer legs. 
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of these effects are found to be of relevance for the motion of a classical light wave around 
the ABCD interferometer loop,[] with explicit calculation (the details will be published 
elsewhere |10|| ) then showing that the ensuing light rays again follow Fig. (1) around the 
interferometer. However, even while there is still a spatial offset A\A 2 = gH 2 /c 2 just as 
before, for light neither a time delay nor any net phase shift is found between the A\C\D 2 
and AB\D 2 paths. However, again as with the neutron case, the spatial offset itself leads to a 
time delay A 2 A/c, so that there is still observable interference. Then, with a replacing mc/h 
in the normalization of the phase shift, we thus find that in traversing the interferometer 
the two beams acquire a final net relative phase shift —aogH 2 /c 2 due entirely to the A 2 A 
segment alone. On recognizing that «o is the value of k° at z = we may set it equal to 
2tt/X where A is the wavelength of the incident beam, to finally obtain for the phase shift 
A<Pcl = —2irgH 2 /\c 2 where CL denotes classical light. Now while H would have to be of 
the order of 10 5 cm for A<Pcl to actually be detectable in a Bragg scattering interferometer 
of the same sensitivity as the COW experiment ,0 nonetheless we can still identify this phase 
shift as an in principle, completely classical effect which reveals the intrinsically global nature 
of classical gravity.0 

Now that we have obtained A^ql = —1 r ngH 2 j\c 2 it is instructive to compare it with 
Acficow- If we introduce the neutron de Broglie wavelength A n = h/mvo we may rewrite 
Acficow in the form —2ngH 2 / \ n VQ. Comparison with A(pcL thus reveals a beautiful example 
of wave particle duality, with the quantum mechanical matter wave inheriting its interference 
aspects from the behavior of the underlying classical wave. Thus even while Aipcow does 
depend on the mass of the neutron]^] its dependence is strictly kinematic with gravity only 

7 As regards the Bragg scattering rules for light in curved space, it is shown in (To) that a wave undergoes 
no change in the magnitude of k° in a Bragg scattering, that it emerges with an angle of reflection equal to 
the angle of incidence, and that the magnitudes of the spatial components of the outgoing momentum take 
whatever values are needed to keep the outgoing wave on the light cone of Eq. (||). Thus, an upward going 
k 3 = «o light wave, for instance, will reflect at C into an initially horizontal light wave with momentum 
k 1 = «o// 1/2 (H). 

8 Whilc this is still sizable for an interferometer, its interest lies in the fact that it allows us to detect, in 
principle at least, the gravitational bending of light using laboratory sized distance scales rather than solar 
system sized ones. Thus it would also be of interest to see what dimension interferometer might serve as a 
gravitational wave detector or be sensitive to any possible neutrino masses in neutrino beam interferometry. 

9 Given our long experience with Newtonian potentials, we are used to the notion that gravity is local, 
with the most distant matter making the smallest contribution at a local point (and none at all if distributed 
spherically), with the non- vanishing of AipcL thus exposing a fundamental difference between Newtonian 
and relativistic gravity. Thus it is of some interest to note that the potential of conformal gravity || takes 
the form V(r) = —MG/r + ~/c 2 r/2 where 7 is a new parameter for sources, to show that there are theories of 
gravity in which even the potentials have global aspects, with the most important gravitational contributions 
at a given point then coming from the most distant objects, a very Machian situation. In fact, it has very 
recently been shown |ll| that because of this global aspect of conformal gravity, the global Hubble flow then 
explicitly modifies the motion of matter within individual galaxies in a way which is then able to account 
for the observed systematics of galactic rotation curves without any apparent need for galactic dark matter. 

10 The reason why A(j)cow actually depends on m at all is that even though the position of the minimum 
of the classical action is independent of m (the equivalence principle) nonetheless the actual value of the 
classical action in this minimum does depend on m (though only as a purely kinematic overall multiplying 
factor for —mc j dr which does not affect the position of the minimum); and even though the value of the 
classical action is not observable classically, nonetheless it is observable quantum mechanically as the phase 
of the wave function whose normalization then explicitly depends (kinematically) on to. 
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coupling via the neutron's de Broglie wavelength, with the COW experiment thus apparently 
being completely compatible with the equivalence principle. 
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Figure Captions 

Figure (1). The paths followed by waves in a COW type interferometer. 
Figure (2). The paths followed by waves in a double slit experiment. 
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Figure 1 




Figure 2 



